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Measurement-based quantum computation utilizes an initial entangled resource state and pro- 
ceeds with subsequent single-qubit measurements. It is implicitly assumed that the interactions 
between qubits can be switched off so that the dynamics of the measured qubits do not affect the 
computation. By proposing a model spin Hamiltonian, we demonstrate that measurement-based 
quantum computation can be achieved on a thermal state with always-on interactions. Moreover, 
computational errors induced by thermal fluctuations can be corrected and thus the computation 
can be executed fault-tolerantly if the temperature is below a threshold value. 
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Introduction.- A quantum computer can solve certain 
problems considered hard for a classical computer with 
an exponential speedup [1]. Standard quantum comput- 
ing uses unitary evolution as a basic mechanism for infor- 
mation processing. Another paradigm is measurement- 
based quantum computing (MBQC), in which one pro- 
cesses quantum information by single-particle opera- 
tions and measurements only, on a nontrivial entangled 
state [2]. Such entangled states serve as universal re- 
sources of MBQC Q. The first identified universal re- 
source was the cluster state. It can be obtained as the 
unique ground state of a Hamiltonian with five-body in- 
teractions ;4] , but can never occur as the unique ground 
state of any two-body Hamiltonian |5j. Fortunately, 
there exist universal resources that are the unique ground 
states of two-body Hamiltonians, albeit with particles of 
local Hilbert space larger than that of a qubit. These 
two-body Hamiltonians include the tricluster model 6 1 , 
an Affleck-Kennedy-Lieb-Tasaki (AKLT)-like model [7[, 
the two dimensional AKLT model [qj9| and a quadratic 
Hamiltonian of continuous variables jlC| . However, in or- 
der to use the ground state of a system as a universal re- 
source, one usually needs to switch off interactions of the 
system sequentially [l,0,&[ll|- Otherwise, the desirable 
quantum correlations could be destroyed due to the time 
evolution of the state via interactions. Therefore, in pre- 
vious proposals, MBQC based on ground states requires 
not only single-particle operations and measurements but 
also a good control of interactions. In this paper, we find 
that it is possible to remove this extra requirement, i.e., 
MBQC can be performed with always-on interactions. 

To this end, we propose a two-dimensional (2D) sys- 
tem and a three-dimensional (3D) system, whose ground 
states are universal resources for MBQC. We show that 
2D and 3D systems can be generalized to a family of sim- 
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FIG. 1: The two-dimensional system composed of spin-3/2 
particles. The system is a hexagonal lattice, where center 
particles (red round) are located on vertices, while bond par- 
ticles (blue ring) are located on edges. There are two kinds of 
interactions between center particles and bond particles, Viinc 
(line) and Vdash (dash line), r denotes the position of a center 
particle, and the vectors between the center particle and its 
three interacting bond particles are 1,2,3 respectively. 



ilar models. These spin models may be realized in physr 



cal systems such as cold atoms [12| , polar molecules 13[ , 
trapped ions [3] and Josephson junction array isj . 
We construct a ground state as a universal resource for 
MBQC by showing that the ground state can be con- 
verted into a cluster state by single-particle operations 
and measurements [^, [l6| . In practice, one obtains a 
thermal state instead of the ground state as a universal 
resource for MBQC. Thus an energy gap is needed to pro- 
tect the state from thermal fluctuations, which is indeed 
the case in our model. However, it is not clear how high 
a temperature can be tolerated before the state would no 
longer be a universal resource of MBQC. Therefore, we 



investigate their thermal states and find that computa- 
tional errors in MBQC induced by thermal fluctuations 
can be corrected as long as the temperature is below a 
certain value. 

2D System.- The 2D system is shown in Fig. [U which 
is a hexagonal lattice with one more particle on each edge. 
The system is composed of spin-3/2 particles, in which 
particles on edges are called bond particles, while others 
are called center particles. Particles are combined by two 
types of interactions 



Mine = AiS^^Al + SyAl + S^A^^), 



(1) 
(2) 



where is the spin operator of the corresponding cen- 
ter particle, and Ag, are operators of the correspond- 
ing bond particle [17]. Operators of bond particles sat- 
isfy commutation relations [I^,J^] = i^ij^aPjI-t, and 
= 3/4, where I,J — A,B, and a,l3,j — x,y,z. 



Therefore, A^ and B^ are two sets of independent spin- 
1/2 operators. 

The Hamiltonian of the system is H = Y^r^r, where 

hr = ■ 'fr, "^r = "^'r-l-l + r+2 + and r 

denotes the position of a center particle, r + a denotes 
the position of one bond particle interacting with the 
center particle r. Here, {a} depends on / e {A, B} as 
shown in Fig. [T] 

Ground state and energy gap.- We can rewrite as 
hr = A(T^2 _ ^2 _ 'f2y2, where I^r = + ^r- Here, 
~§ r, ~^ r and T^r all satisfy commutation relations of spin 
operators. Therefore, 



+ l)-5r(5r + l)-/r(/r + l)], (3) 



where Sr = 3/2 and /r = 1/2 or 3/2. When 4 = 1/2, 
Tr = 1, 2. When 4 = 3/2, = 0, 1, 2, 3. One can get the 
minimum energy by taking Ij. — Z/2 and Tr — 0, which 
means the ground state, Ig)^., of hj. has a total spin 0. 
The energy difference between the ground state and the 
first excited state is A. Because these are independent 
with each other, the ground state of the whole system is 
|G) = {^j. |g)j. and protected by an energy gap A. The 
energy gap only depends on the interaction constant, and 
does not vanish in the thermodynamic limit. 

POVM and GHZ state - As the first step of MBQC on 
the ground state, the POVM I = Y.a=x,v,z F°''' F°' is per- 
formed on center particles. Here, F" = {S"^ — 1/4)/ \/6, 
which projects the center spin into the subspace spanned 
by two states with maximum spin component in the a di- 
rection. Because the ground state \g)^ has a total spin 0, 
all three spin-/r+a are antiparallel with the center spin- 
Sr- Therefore, the POVM projects the state |g)j. into a 
GHZ state, e.g., for the outcome z, the output state is 
|ghz)^ = (10000) + |Tlll))/v^, where |0) = -\S^ = 3/2), 
|1) = \S^ = —3/2) are the state of the center spin, and 
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FIG. 2: The three-dimensional system, (a) The elementary 
cell of the system. The system is composed of spin-2 parti- 
cles and spin-3/2 particles. Spin-2 particles are center par- 
ticles (red round), and spin-3/2 particles are bond particles 
(blue ring), (b) The elementary cell of the three dimensional 
topology-protected cluster state, which can be prepared by 
single-particle operations and measurements on the ground 
state of the 3D system, (c) Three directions for fc < 3 of the 
POVM on spin-2 particles, which are orthogonal with each 
other and passing face centers of a cube, (d) Four directions 
for > 4 of the POVM on spin-2 particles, which are along 
body diagonals of the same cube. 



|0) (|1)) is the eigenstate of I^+a with eigenvalue —1/2 
(1/2). The state \g)^ is an isotropic state. Therefore, all 
outcomes are equivalent to the outcome z up to a set of 
single-particle operations U{a) = exp[i (a)], where 
a is the outcome and lt{a) = axzarcsin(|Sxz|)/|Sxz|. 
Then, the state of the whole system after POVMs and 
single-particle operations is |{ghz}) = Ighz)^., which 
can also be described by a set of stabilizers, Wr — 

Vr na=i.2.3 2^r+a and lFr,r+a = 2^r-?^^+a fo^' aU r and 
a. |{ghz}) is the eigenstate with eigenvalue 1 of all of 
these stabilizers. Here, X, Y and Z are Pauli operators 
of the qubit {|0), |T)}. 

Cluster state and universality of the ground state - By 
measuring physical quantities A^B^ and A^Bl on bond 
particles, the state |{ghz}) can be projected 18, 1^ (or 
"fused") into a hexagonal cluster state, which has the 
same lattice with center particles. Eigenstates of A^B^ 
and A^B^, which are measurement basis, can be found 
in Ref. [Uj . Considering a product of stabilizers 



Wr 



n 



r+2a,r+a 



Zr+2a(4Ar+a-Sr+a)) 



(4) 

one can get a new stabilizer by replacing A^B^ with out- 
comes. In Eq. (HI), we have taken the case /^+a = ^r-i-a as 
an example, and the result is the same for 7^+a = ^r+a- 
The new stabilizers define a hexagonal cluster state on 
center particles up to a Pauli frame, which can be cor- 
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rected by single-particle operations The hexagonal 
cluster state is a universal resource for MBQC [3|- Then, 
universal MBQC can be performed on center particles. 

3D system and topology-protected cluster state - Fol- 
lowing the idea of 2D system, we propose a 3D sys- 
tem, whose ground state is also a universal resource 
for MBQC. The system is shown in Fig. [2] (a), which 
is composed by spin-2 particles and spin-3/2 particles, 
where center particles are spin-2 particles and bond par- 
ticles are spin-3/2 particles. The interactions between 
particles are the same as Eq. (jlj and ([2]). Therefore, 
the Hamiltonian of the 3D system has the same form 
as 2D system, H = ^^.h^, = ^ • where 

7"r = ~fr+i + ~fr+2 + ^ r+3 + ~t v+i- Here, {r + a} 
denote four bond particles interacting with the center 
particle r. 

In the 3D system, one can get the minimum energy of 
hr by taking 1^ — 2 and Tr = in Eq. ([3]). Therefore, in 
the 3D system, the ground state of each is an isotropic 
state with a total spin 0. The energy difference between 
the ground state and the first excited state is A, which 
means the 3D system is also gapped. 

The ground state of the 3D system can be reduced 
to a 3D cluster state, as shown in Fig. [5] (b). Firstly, 
center particles are measured as I = J2k=i P'' {^k)F{ak), 
which is a POVM with seven outcomes. Here, F{ak) — 
V7VfcP(Sfc), and 



P(S) = |S;2) (S;2H-|a;-2) (a;-2| 



(5) 



projects the center spin into the subspace spanned by 
two states with maximum spin component in the a di- 
rection. \a;m) is the eigenstate of a ■ ~§ c with eigen- 
value m. Nk = 1/3 for fc < 3 and Nk = 3/8 for 
fc > 4. The seven directions are shown in Fig. [2] (c) 
and (d). Because four spins {7r-i-a} are all antiparallel 
with the center spin-S'r, the output states of the POVM 
are GHZ states. These GHZ states are equivalent to the 
GHZ state of outcome z, up to single-particle operations 
U{a). Therefore, POVMs on center particles, with U{a) 
together, can transform the ground state to a state sta- 

bihzed by Wr = Xr S^r+a and Wr,r+a - 2ZJ^+^, 

where a = 1, 2, 3, 4. 

Measuring physical quantities A^B^ and A^B^, one 
can generate a new set of stabilizers Ha ^r+2a, which 
defines a 3D cluster state on center particles, as shown 
in Fig. [2] (b). On the 3D cluster state, quantum cor- 
relations are protected topologically and fault tolerant 
quantum computing can be simulated using topological 
error correction 20]. 

Thermal computational errors and error correction - 
We have proved the ground states of 2D and 3D systems 
are universal resources for MBQC. However, in practice, 
a system cannot reach the exact ground state, but rather 
a thermal state at finite temperature. Thermal fluctu- 
ations can reduce the quantum correlations on ground 
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FIG. 3: Error probabilities e on a GHZ state and p on the 
cluster state as functions of the temperature. For 2D system 
the GHZ state is a four-qubit state, and for 3D system the 
GHZ state is a five-qubit state, e are almost the same for 
2D and 3D systems when the temperature T < A. When 
T/A = 0.2, p = 3%, which means errors induced by a lower 
temperature are tolerable by using the topological error cor- 
rection algorithm on the 3D system. 



states and induce computational errors on the cluster 
state, which will be used for MBQC. The thermal state 
is the Gibbs state p = Z~^e~^" , where Z = tre"^^, 
/3 = l/T is the temperature, and p can be rewritten as 
p — pr- Here pr = Z~^e~^''"' is the Gibbs state of hr- 
After the POVM and U{a), the state pr is transformed 
into CTr = FprFytr{FprF^), where F ^ F"^ for the 2D 
system and F = F{z) for the 3D system. At an absolute 
zero temperature, a-r — |ghz) (ghzj^. is the desired GHZ 
state. Here, Ighz)^. is a GHZ state of four qubits for the 
2D system and five qubits for the 3D system. 

The post-POVM state CTr at finite T is only approxi- 
mately a GHZ state, i.e., is equivalent to a perfect GHZ 
state affected by errors. The probability of an error oc- 
curring on the post-POVM state is e = 1 — F, where 
F ~ tr(f7r Ighz) (ghzlr) is the fidelity of the GHZ state, 
as shown in Fig. |31 Those errors are propagated under 
the measurements of the bond particles and subsequent 
correction operations 2l|. The resulting error superop- 
erators act on the yet unmeasured center particles, and 
have the following properties: (I) there is one indepen- 
dent error superoperator Sr for every r, (2) £r acts at the 
locations r and {r -I- 2a, Va}. Where the center particles 
are measured in the X-basis for the purpose of topologi- 
cal error correction on the 3D cluster state (in most of the 
cluster), there arise two further simplifications: (3) All 
errors are equivalent to Z-errors or the identity, and (4) 
Correlations between errors on neighboring center parti- 
cles can be discarded. The latter arises because errors at 
r and at r -I- 2a are corrected by different error-correction 
procedures running independently of another [20| . 
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On the 3D cluster state, for each r, the resulting error 
is £r = Eio E2, with El = {I — pi) + pi[Zr] and E2 = 

(1 -P2 -Pz) +P2/4 Ea[^r+2a] +P3/6 Ea,a' [^r+2a^r+2a' ] ■ 

Therein, the error probabilities pi, P2 and p^ depend 
on the temperature T. If ^ Pi,P2 ^ 1, then 
the local errors are almost independent and the error 
level is described by an effective local error probabil- 
ity p ~ pi + P2 + 2p3. Error-correction is possible if 
p < 3% 221, which translates into a threshold tempera- 
ture Tt = 0.2A; See Fig. H At that point, pi,p2 ~ 10"^ 
and P3 = 10^^, justifying the assumption of uncorrelated 
local errors. 

MBQC with always-on interactions - In practical appli- 
cation, one can convert the initial state, usually a thermal 
state, to a cluster state one qubit at a time. Once we need 
the qubit r, we can apply POVMs on the center particle 
r and its neighboring center particles {r-|-2a}. Based on 
outcomes of POVMs, single-particle operations U (a) are 
chosen. Then, bond particles {r -I- a} are measured, and 
outcomes are used to correct the Pauli frame of qubit r. 
No further operation is needed on any other particle in 
order to convert the center particle r to a qubit on the 
cluster state. 

With always-on interactions, we need to consider the 
time evolution driven by the time-independent Hamil- 
tonian. Since the initial state is not converted to the 
cluster state simultaneously, there are some untouched 
particles. They remain in the initial thermal state, which 
is close to the ground state due to the existence of the en- 
ergy gap. Other particles evolve with the Hamiltonian, 
and their quantum correlations will be changed under 
time evolution. Fortunately, the time evolution is peri- 
odic with a period 47r/A for the 2D system and 2tt/A 
for the 3D system. Then, one can perform operations on 
these particles at the revival time of quantum correla- 
tions, t = 4n7r/A and t = 2r7,7r/A for 2D and 3D systems 
respectively, where t = is the time of the first oper- 
ation on the particle and its interaction particles, and 
n = 0, 1, 2, . . .. If we assume only one operation can be 
performed on each particle at one revival time, particles 
can be measured out before n = 6. Therefore, the MBQC 
can be performed on our proposed systems with always- 
on interactions. Here, operations are required to occur 
precisely. We remark that errors due to timing impreci- 
sion can also be analyzed similarly to thermal errors. 

Discussion - In summary, we proposed a 2D and a 3D 
gapped system, whose ground state is entangled based 
on a factorized Hamiltonian. With a factorized Hamil- 
tonian, quantum computing can be performed without 
the need to switch off interactions. The ground state can 
be reduced to a deterministic cluster state, in contrast to 
AKLT like systems where cluster states are obtained with 
stochastic structures. Errors induced by thermal fluctu- 
ations can be corrected as long as the temperature is 
below a critical threshold. There are other choices of 
and that satisfy conditions of spin-1/2 operators. By 



replacing the center particle with different spin systems, 
i.e. spin-m/2 particles, one can get different spatial con- 
nectivities that each qubit is connected to m other qubits 
in the cluster state. Thereby, it can be generalized to 3D 
and more complicated configurations. 
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